NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


WARTIME  REPORT 


ORIGINALLY  ISSUED 
May  1944  as 

Restricted  Bulletin  L4E29 


HEAT -CAPACITY  LAG  IN  TURBINE  -WGEKOSG  FLUIDS 
By  Arthur  Kantrowltz  and  Paul  ¥.  Huber 


Langley  Memorial  Aeronautical  Laboratory 
Langley  Field,  Va. 


WASHINGTON 


^SEaF***** 

Ur Ztey  Field,  V*. 


NACA  WARTIME  REPORTS  are  reprints  of  papers  originally  issued  to  provide  rapid  distribution  of 
advance  research  results  to  an  authorized  group  requiring  them  for  the  war  effort.  They  were  pre¬ 
viously  held  under  a  security  status  but  are  now  unclassified.  Some  of  these  reports  were  not  tech¬ 
nically  edited.  All  have  been  reproduced  without  change  in  order  to  expedite  general  distribution. 


3  1176  01364  8549 

NACA  RB  No.  ll|E29 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


RESTRICTED  BULLETIN 

HEAT-CAPACITY  LAC  IN  TURBINE-WORKING  FLUIDS 


By  Arthur  Kantrowitz  and  Paul  W.  Huber 


SUMMARY 


Gases  require  a  finite  time  -  the  relaxation  time  - 
to  adjust  the  vibrational  part  of  their  internal  energy 
to  a  change  in  temperature.  If  changes  In  temperature 
occur  in  gas  flow  In  times  of  the  order  of  or  shorter  than 
the  relaxation  time,  entropy  Increases  occur.  In  order 
to  evaluate  partially  the  importance  of  these  heat- 
capacity  lags  In  turbines,  the  relaxation  times  of  steam 
and  nitrogen  have  been  studied  by  the  total-head-defeot 
method. 

The  relaxation  time  of  steam  was  found  to  be  of  the 

q 

order  of  3  *  10“°  second.  It  would  be  expected  from 
this  result  that  the  vibrational  heat  capacity  of  steam 
would  adjust  to  moat  of  the  temperature  changes  In  en¬ 
gineering  practice  rapidly  enough  that  little  dissipation 
would  be  anticipated.  Water  vapor  Is  known  to  be  an 
efficient  catalyst  for  combustion  processes  and  the  rapid 
adjustment  of  vibrational  heat  capacity  may  be  Important 
In  its  catalytic  actions. 

Adjustment  of  the  vibrational  heat  content  of 
nitrogen  to  a  change  in  temperature  was  found  to  be  a 
slow  process  and  principally  dependent  on  the  action  of 
impurities.  Measurements  were  made  of  mixtures  of  steam 
with  nitrogen  and  it  was  found  that  about  1  collision 
in  30,000  between  steam  and  nitrogen  molecules  is  effec¬ 
tive  in  adjusting  the  vibrational  heat  content  of  nitrogen. 
Nitrogen  at  1100°  F  absolute  with  3  percent  water  vapor 

thus  would  have  a  relaxation  time  of  1.5  x  10”^"  second  at 
atmospheric  pressure  so  that  heat-capacity  lag  in  nitrogen 
should  be  considered  in  the  design  of  gas  turbines. 


INTRODUCTION 


Turbine -working  fluids  such  as  steam,  air,  and  ex¬ 
haust  gas  have  appreciable  vibrational  heat  capacity  at 
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high  temperatures.  If  these  gases  have  relaxation  times 
of  the  order  of  or  shorter  than  the  intervals  during 
which  temperature  changes  occur  in  the  gas,  losses  attrib¬ 
utable  to  heat-capacity  lag  must  be  anticipated.  Uhfortu 
nately  no  measurements  of  the  relaxation  times  of  the 
most  important  turbine -working  fluids  have  previously  been 
made.  Measurements  of  the  relaxation  time  of  steam  and 
nitrogen  are  reported  in  the  present  paper. 

The  measurements  were  made  by  the  total-head-defect 
method,  which  was  introduced  in  reference  1  and  has  been 
herein  extended  to  permit  measurement  in  gases  with  long 
relaxation  times. 

The  symbols  used  throughout  the  present  paper  are 
defined  in  appendix  A. 


TEST  METHODS 


A  drawing  of  the  basic  apparatus  used  to  obtain  the 
measurements  presented  herein  is  shown  in  figure  9  of 
reference  1  and  the  apparatus  Is  described  in  detail  on 
page  21  In  reference  1.  An  outer  lining  on  the  chamber 
and  fine-mesh  screens  Inside  the  chamber  were  incorporated 
in  the  apparatus  used  In  these  experiments.  An  electric 
boiler  and  superheater  were  constructed  which  produced 
enough  steam  for  an  hour's  run.  A  heated  manometer 
system  was  used  in  the  steam  measurements  to  avoid  con¬ 
densation  difficulties.  The  boiler  was  also  used  to  add 
water  vapor  to  the  nitrogen  and  a  dew-point  instrument  was 
used  during  the  nitrogen  tests  to  measure  the  water-vapor 
content. 

The  steam  boiler  was  scrubbed  and  cleaned  a  number 
of  times  and  then  heated  In  a  vacuum  for  some  time  before 
it  was  used.  Distilled  water  was  used  at  all  times  In 
the  boiler  and  only  the  middle  two-thirds  of  the  steam 
generated  was  used  in  these  experiments.  Oases  Insoluble 
In  water  were  found  In  the  steam  to  about  one  part  in  a 
million. 

The  nitrogen  used  was  commercial  dry,  99  percent  pure, 
obtained  in  high-pressure  cylinders.  Analysis  of  a 
sample  of  this  nitrogen  showed  that  it  contained  less  than 
0.1  percent  oxygen  and  less  than  0.25  percent  carbon 
dioxide.  Passage  of  the  gas  through  the  dew-point  instru¬ 
ment  indicated  that  roughly  0.05  percent  water  vapor  was 
present. 
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HEAT-CAPACITY  LAO  IN  STEAM 


Total -head  defects  of  steam  were  measured  at  three 
different  temperatures  and  with  two  Impact  tubes  over  a 
range  of  chamber  pressure':  The  total-head  defects 

measured  are  plotted  against  chamber  pressure  In  figures  1 
and  2.  The  heated  total-head-defect  nenometer  used  made 
It  difficult  to  obtain  an  accurate  zero  reading.  The 
zero  readings  were  therefore  obtained  by  extrapolating 
the  experimental  curve  to  zero  chamber  pressure.  The 
results  obtained  are  compared  with  the  theory  of  refer¬ 
ence  1  in  figures  1  and  2.  It  will  be  seen  from  fig¬ 
ure  1  that  the  total-head  defect  is  roughly  inversely 
proportional  to  the  impact-tube  diameter.  This  fact  and 
the  fact  that  the  defects  observed  were  small  means  that 
steam  has  a  relaxation  time  which  is  short  compared  with 
the  compression  time  at  the  nose  of  the  impact  tubes. 

The  relaxation  times  were  then  calculated  by  employing 
the  theory  of  total-head  defects  due  to  compressions  at 
the  nose  of  an  impact  tube  (reference  1).  By  choosing 
three  typical  experimental  points,  the  relaxation  times 
shown  in  figure  3  were  found.  By  using  these  experi¬ 
mental  values  of  relaxation  time  and  the  theory  of  refer¬ 
ence  1  and  presuming  that  the  collisions  responsible  for 
the  adjustment  of  vibrational  energy  are  bimolecular 
(collisions  involving  two  molecules;  relaxation  time  t 
inversely  proportional  to  pressure),  the  theoretical  total- 
head  defects  shown  in  figures  1  and  2  were  computed.  The 
agreement  of  the  theoretical  variations  with  chamber  pres¬ 
sure  and  impact-tube  diameter  and  these  measurements  may 
be  taken  as  an  index  of  the  accuracy  of  these  results, 
which  is  probably  of  the  order  of  ±20  percent  for  relaxa¬ 
tion  time.  Inasmuch  as  the  steam  employed  probably  had 
only  slight  impurities  and  the  adjustment  involved  only  a 
few  hundred  collisions  per  molecule,  these  measurements 
maybe  presumed  representative  of  pure  steam.  Some 
experimental  relaxation  times  for  steam  are  shown  in 
table  I.  It  can  be  seen  that  a  very  short  time  -  of  the 

-8 

order  of  3  x  10  second  -  is  required  for  adjustment. 

The  catalytic  action  of  steam  may  be  partly  attributed  to 
this  rapid  interchange  of  energy. 

The  results  shown  in  figures  1  and  2  are  next  com¬ 
pared  with  the  theory  of  Landau  and  Teller  (reference  2). 
This  theory  indicates  that  the  logarithm  of  the  number  of 
collisions  necessary  to  transmit  a  vibrational  quantum  is 
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inversely  proportional  to  the  cube  root  of  the  absolute 
temperature  T.  The  number  of  collisions  required  for 

_1 

the  relaxation  of  steam  is  plotted  against  T  ^  in 
figure  ij..  The  results  presented  herein  as  well  as  those 
of  other  recent  investigations  (see  reference  1)  have 
been  shown  to  fit  the  theory  of  reference  2  within  experi¬ 
mental  accuracy. 


HEAT-CAPACITY  LAG  IN  NITROGEN 


Total-head  defects  of  nitrogen  were  measured  at 
different  temperatures  over  a  range  of  chamber  pressure  1 
with  essentially  the  same  apparatus  as  for  the  measure¬ 
ments  with  steam.  No  change  in  total-head  defect  was 
found  with  change  of  diameter  of  the  Impact  tube.  This 
result  Indicated  that  the  relaxation  time  of  nitrogen 
was  long  compared  with  the  compression  time  at  the  nose 
of  the  impact  tubes  usefl.  The  addition  of  small  amounts 
of  Impurities  in  order  to  shorten  the  .relaxation  time 
resulted  In  much  larger  total-head  defects  and  further 
established  the  fact  that  the  relaxation  time  of  nitrogen 
Is  of  the  same  order  as  the  time  required  for  the  expansion 
through  the  nozzle.  The  total-head  defects  were  there¬ 
fore  computed  for  a  flow  of  this  type  by  extending  the 
theory  of  reference  1.  This  calculation  is  reported  In 
appendix  B  and  was  used  to  find  the  relaxation  time  of 
nitrogen. 

The  relaxation  times  measured  for  commercial  nitrogen 
and  nitrogen  with  various  amounts  of  added  water  vapor 
are  shown  in  table  I.  It  can  be  seen  that,  for  the 
relatively  pure  gas,  a  time  longer  than  0.005  second  Is 
required  whereas.  In  a  steam  and  nitrogen  mixture, 
about  1  collision  in  30,000  between  steam  and  nitrogen 
molecules  is  effective.  Nitrogen  at  1100°  P  absolute 
with  3  percent  water  vapor  therefore  would  require 

about  1.5  x  10"^-  second  at  atmospheric  pressure  to  adjust 
its  vibrational  heat  capacity.  It  can  also  be  seen  that 
the  addition  of  water  vapor  shortens  the  relaxation  time 
by  a  factor  roughly  proportional  to  the  amount  added. 

As  has  been  observed  with  some  other  gases  (reference  3)* 
the  dependence  of  the  number  of  collisions  for  relaxation 
on  temperature  is  much  smaller  for  a  gas  with  Impurities 
than  for  a  relatively  pure  gas. 
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A  short  investigation  was  made  to  determine  whether 
carbon  dioxide  in  the  presence  of  water  vapor  is  an  ef¬ 
fective  catalyst  for  the  transformation  of  the  vibra¬ 
tional  energy  of  nitrogen  into  translational  energy.  The 
speculation  was  introduced  by  Mr.  L.  R.  Turner  of  the 
Aircraft  Engine  Research  Laboratory  and  was  based  on  the 
fact  that  the  frequency  of  one  of  the  vibrational  modes 
of  CC>2  and  that  of  nitrogen  are  about  in  resonance. 

Total -head  defects  were  obtained  with  a  1. 5 -inch  nozzle 
and  0,060-inch  impact  tube  on  a  mixture  of  nitrogen  with 
3-percent  water-vapor  content  at  a  temperature  of  about 
800°  P.  Amounts  of  CO2  up  to  5  percent  were  added  to 
this  mixture  and  only  small  changes  (less  than  one- 
twentieth  of  the  defect)  were  noted.  In  a  mixture  of 
this  kind,  CO2  is  therefore  considered  to  have  only  a 

small  effect  on  the  relaxation  time  of  nitrogen. 


APPLICATIONS 


The  long  relaxation  time  found  for  nitrogen  means 
that  the  vibrational  energy  of  this  gas  lags  behind  the 
rapid  temperature  changes  which  occur  in  many  engineering 
applications.  As  an  illustration  of  the  type  of  effect 
to  be  anticipated,  the  expansion  of  a  mixture  of  hot  gases 
through  a  converging-diverging  nozzle  from  a  pressure  of 
I4.  atmospheres  to  a  pressure  of  0.5  atmosphere  has  been 
analyzed.  The  mixture  Is  composed  of  JO  percent  air  and 
30  percent  exhaust  gas  and  components  other  than  water 
vapor  are  assumed  to  have  no  effect  on  the  relaxation 
time  of  nitrogen.  (CO2  has  been  shown  to  have  little 
effect.)  The  relaxation  times  of  the  components  other 
than  nitrogen  are  assumed  to  be  zero.  In  order  to  sim¬ 
plify  the  calculations,  the  nozzle  was  assumed  to  be  so 
designed  that  the  translational  temperature  decreased 
linearly  with  time.  The  translational  and  vibrational 
temperatures  found  are  plotted  against  time  In  figure  5* 
(Translational  temperature  is  defined  as  the  temperature 
that  would  be  measured  In  the  gas  at  equilibrium  having 
the  same  amount  of  translational  kinetic  energy.  Vibra¬ 
tional  temperature  Is  defined  as  the  temperature  that 
would  be  measured  in  the  gas  at  equilibrium  having  the 
same  amount  of  vibrational  energy.)  It  is  seen  that  very 
little  adjustment  of  the  vibrational  heat  capacity  occurs 
in  the  nozzle.  The  nozzle  is  0.1  foot  long,  its  area 
ratio  in  the  diverging  part  is  designed  for  an  8:1  expansion 
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ratio  when  lsentropic  flow  Is  presumed,  and  the  Initial 
temperature  is  2lj.60o  F  absolute. 

The  results  of  approximate  calculations  are  presented 
in  the  following  table: 


Losses  in 
mechanical 
energy 
available 
(percent) 

Pres¬ 

sure 

ratio 

across 

nozzle 

Mass  flow 
(s]ugs/&ec/feq  ft) 

Velocity 

at 

nozzle 

exit 

(fps) 

lsentropic 

flow 

0 

8 

2.76 

5723 

Plow  Including 
IjLeat-capac  ity 

2 

8.i|0 

2.79 

3712 

CONCLUSIONS 


The  vibrational  heat  capacity  of  steam  was  found  to 
be  adjusted  very  rapidly  by  bimolecular  collisions  in  the 
gas.  A  few  hundred  collisions  per  molecule  or  a  time  of 
the  order  of  3  x  10”“  second  Is  required.  The  rapidity 
of  this  adjustment  means  that  the  vibrational  heat  capa¬ 
city  of  steam  follows  the  temperature  changes  in  most 
flow  processes.  This  rapid  interchange  may  be  important 
in  the  catalytic  actions  of  steam. 

The  vibrational  heat  capacity  of  nitrogen  was  found 
to  be  adjusted  only  slowly  in  the  pure  gas,  a  time  longer 
than  0.005  second  being  required.  The  adjustment  is 
materially  accelerated  by  the  presence  of  water  vapor, 
about  1  collision  In  JO ,000  between  steam  and  nitrogen 
molecules  being  effective.  Nitrogen  at  1100°  F  absolute 
and  atmospheric  pressure  with  a  5-percent  water-vapor 
content  would  have  a  relaxation  time  of  1.5  x  10 
second.  The  water  vapor  in  exhaust  gas  thus  does  not 
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accelerate  the  adjustment  of  nitrogen  enough  that  its 
vibrational  heat  capacity  can  follow  the  rapid  tempera 
ture  changes  In  the  flow  through  the  gas  turbines. 


Langley  Memorial  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va. 
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APPENDIX  A 

SYMBOLS 


nozzle  length 
diameter  of  Impact  tube 
final  velocity  in  nozzle 
velocity 

nondime ns lonal  velocity  (u/U) 
time 


nondimens ional  time 


(i/u  or  d/u) 


relaxation  time  at  average  temperature  and  pres¬ 
sure  of  flow 

nondimens ional  constant  of  flow  referred  to  impact 
tube  (d/u-r) 

nondimensional  constant  of  flow  referred  to 
nozzle  (I/Ut) 

excess  of  internal  energy  over  that  of  equilibrium 
at  the  translational  temperature  (nondimensional) 

nondimensional  entropy  increase 

radius  of  nozzle  cross  section 

distance  along  nozzle  axis 

temperature 
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APPENDIX  B 

MEASUREMENTS  OP  HEAT - CAP AC  IT Y  LAC  FOR  GASES  WITH  LONG 
RELAXATION  TIMES  BY  MEANS  OP  A  NOZZLE  OP  SPECIAL  DESIGN 


In  the  measurements  of  heat-capacity  lag  of  gases 
with  the  apparatus  shown  in  figure  9  of  reference  1,  a 
nozzle  is  employed  in  which  the  gases  expand  and  acceler¬ 
ate  before  meeting  the  impact  tube.  If  the  relaxation 
time  of  a  gas  is  long  enough  to  be  of  the  same  order  as 
the  expansion  time  through  this  nozzle,  an  entropy  in¬ 
crease  occurs  in  the  nozzle .  For  ease  of  calculation, 
the  nozzle  used  in  these  tests  was  so  designed  that 


-  =  Constant  =  X 

dt 


(1) 


The  shape  of  a  nozzle,  of  circular  cross  section, 
to  meet  the  condition  of  equation  (1)  can  be  readily 
found  if  low- velocity  flow  -  that  is,  constant  density  - 
is  assumed.  When  the  velocity  is  assumed  constant 
across  the  cross  section,  continuity  yields 

ur^  =  Constant  (2) 

Also, 

u  dt  =  dx  (5) 

Combining  equations  (1)  to  (3)  gives 

r^x  =  Constant  (lj.) 

Equation  (lj.)  must  be  corrected  for  the  presence  of  a 
boundary  layer  along  the  walls  of  the  nozzle.  The  dis¬ 
placement  thickness  of  this  boundary  layer  -  that  is,  the 
distance  the  streamlines  just  outside  the  boundary  layer 
are  displaced  by  its  presence  -  can  be  found  by  the 
methods  of  reference  I}..  It  was  assumed  that  this  boundary 
layer  would  be  laminar.  Equation  (1)  of  reference  I4.  was 
thus  used  to  calculate  the  displacement  thickness.  By 
assuming  typical  flow  conditions,  this  distance  was  then 
added  to  the  nozzle  radius.  The  radius  of  the  nozzle 
cross  section  is  plotted  against  the  dlstanoe  along  the 
center  line  In  figure  10  of  reference  1. 
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The  entropy  Increase  In  a  nozzle  of  this  design  is 
now  calculated.  It  can  be  shown  from  definition  that 
X  =  2/3.  Prom  equation  (11)  of  reference  1, 

€  »  20 Ktft*  dt,  + 


Constant^ 


where  t'  =  0  is  taken  at  the  start  of  the  expansion 
when  c*  =  0.  The  constant  is  then  zero  and  integration 
yields 


(5) 


From  equation  (16)  of  reference  1  for  the  change  in 
entropy  as  t'  changes  from  0  to  t*max 


Prom  equation  (1),  the  nondimensional  time  t’max  for 
the  gas  to  pass  through  the  nozzle  is  ^  =  \\  hence. 


il  16 
AS.N  =  -y= - : 


5%  9%e 


L 


This  entropy  increase  in  the  nozzle  must  be  added  to 
the  increase  at  the  impact  tube  in  order  to  obtain  the 
total  entropy  increase  in  an  expansion  through  a  nozzle 
of  this  type  followed  by  a  compression  at  the  nose  of  a 
source-shape  impact  tube. 

For  the  calculations  of  the  entropy  increase  at  the 
impact  tube  in  reference  1,  it  was  presumed  that  c*  =  0 
initially.  This  assumption  does  not  hold,  however,  for 
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a  gaa  with  a  long  relaxation  time  whloh  has  reached  a 
tianequi  librium  oondltlon  In  the  nozzle;  consequently,  the 
equations  of  reference  1  must  be  modified  to  fit  this 
case  as  follows;  In  equation  (22)  of  reference  1,  for 
the  variation  of  e*  with  time  at  the  impact  tube,  the 
constant  was  zero  because  e*  was  assumed  zero  before 
compression.  For  a  gas  with  a  long  relaxation  time, 
the  constant  becomes  equation  (5)  with  t’  =  3/2,  which 
is  the  value  of  e»  at  the  end  of  expansion.  This 
assumes  that  the  impact  tube  is  placed  so  close  to  the 
end  of  the  nozzle  that  compression  starts  immediately 
after  the  gas  leaves  the  nozzle,  Equation  (22)  of 
reference  1  now  becomes 


c'(t') 


l6u»2(l  -  u» )2  eKt»dt' 


By  integrating  equation  (6)  by  Simpson’s  rule  and 
equation  (16)  of  reference  1  by  the  planimeter,  both 
between  t’  =  0  and  t’  =  «,  the  entropy  increase  in  the 
compression  at  the  impact  tube  is  obtained  for  a  gas  with 
any  relaxation  time.  The  total  entropy  increase  becomes 

AS’t  =  As 1  tube  +  as’n 

which  is  plotted  against  K  in  figure  6  and  compared  with 
the  theory  of  reference  1, 
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TABLE  I.-  EXPERIMENTAL  RELAXATION  TIMES  OP  STEAM  AND  NITROGEN 


Typ9  and  per- 

Diam.  of 
Impact  tube 
(in.) 

Relaxation 

Molecular 

Type  of  colli  ai  on 

cent  of 

■Temperature 

time  at 

collisions  re- 

considered 

impurities 

(°P  abs. ) 

1  atm,  t 

quired  for 

effective  for 

in  gas 

(sec) 

adjustment 

adjustment 

Steam 

Negligible 

0.017 

wmm 

3.7  x  io“8 

4  x  102 

Bimolecular 

Do—— 

.017 

2.6 

2.6 

Do. 

Do - 

.030 

3.1 

3.1 

Do. 

Do—— 

.017 

2.1 

1.9 

Do. 

Nitrogen 

0.05  e2o 

0.060 

5.9  x  io-3 

a3.1  x  io7 

Bimolecular 

.05 

.060 

12  2k 

2.9 

a1.4 

Do. 

.05 

.060 

1236 

2.7 

al-3 

Do. 

.05 

.060 

1250 

3.4 

al.6 

Do. 

.05 

.060 

1370 

2.6  . 

ai  2 

3  To  x  io4 

Do. 

.7 

.060 

6.2  x  10-4 

N2  with  HgO 

1.0 

.060 

1281 

3.9 

2.5 

Do. 

.060 

3.5 

3.5 

Do. 

l.k 

.060 

1277 

2.8 

2.7 

Do. 

1.4 

.017 

1319 

2.6 

2.4 

Do. 

2.0 

.017 

1034 

2 .3 

3.4 

Do. 

3-0 

3.5 

.017 

.017 

1.4 

1.2 

2.9 

3.6 

Do. 

Do. 

Adjustment  computed  as  due  to  collisions  of  nitrogen  with  nitrogen  even  though  It  .  h 
Is  thought  that  small  amount  of  water  vapor  present  was  important.  01 
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Figure  5.-  Theoretical  departure  of  vibrational  temperature  of  nitrogen  from  trandlational 
temperature  in  the  flow  of  exhaust  gas  through  a  1.2-lnch  nozzle  having  an  8:1  expansion  w 
ratio. 
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Figure  6  -  The  total  nonchmensional  entropy  increase  m  nozz/e  ana 


impact  tube.  AS'r  =  A  4'^  +  A  ;  K  K„  =■  94.1  K-^K. 
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